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Abstract 

We employ variational Monte Carlo methods to study the transition to 
strange matter in a simple one-dimensional string-flip model with two fla- 
vors and two colors of quarks. The dynamics of the system are described in 
terms of a many-body potential that confines quarks within hadrons, yet en- 
ables the hadrons to separate without generating unphysical long-range van 
der Waals forces. The model has "natural" low- and high-density limits: it 
behaves as a system of isolated hadrons at low density and as a Fermi gas 
of quarks at high density. We show that the system exhibits a transition to 
strange matter characterized by an increase in the length-scale for confine- 
ment. Yet the small increase at the transition region — of only ten percent 

- suggests that clustering correlations remain strong well into the strange- 
matter domain. Our results put into question descriptions of strange matter 
in terms of noninteracting, or weakly interacting, quarks. 

PACS number(s): 12.39.-x,24.85.+p,26.60.+c 
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I. INTRODUCTION 



The possible existence of stable strange matter — more stable even than 56 Fe — is 
fascinating indeed That this scenario is plausible at high baryon density is simple 
to understand: for two-flavor (up and down) matter the Fermi energy becomes so large 
that the addition of strange quarks becomes energetically favorable, in spite of the larger 
strange-quark mass. The existence of strange matter should have a dramatic impact on fields 
ranging from astrophysics and cosmology to nuclear and particle physics. Indeed, it is widely 
accepted that the high-density environment at the core of neutron stars constitutes a fruitful 
ground for the formation of quark matter. Moreover, the recent proposal of a "spin-up" phase 
in the timing structure of pulsars, speculated to signal the transition from hadronic to quark 
matter, makes the detection of this exotic phase of matter within observational reach [F2J. 
These astronomical observations complement terrestrial searches for strange matter which 
have improved considerably, and will continue to do so, with the commission of new powerful 
relativistic heavy- ion colliders |J. 

The theoretical study of strange matter has traditionally been conducted in two differ- 
ent pictures: one using a hadronic model [|J - similar to ordinary nuclei — where the 
fundamental degrees of freedom are mesons and baryons, another using a quark model |5]|| 
consisting of massless, noninteracting quarks confined inside a bag (note that the inclusion 
of a finite strange-quark mass and a perturbative treatment of QCD does not seem to change 
the qualitative picture). Presumably a description of strange matter in terms of mesons and 
baryons is well motivated in the low-density regime where clustering correlations remain 
important. Alternatively, strange matter viewed in terms of a Fermi gas of quarks might be 
appropriate in the high-density domain, as now the average interquark separation becomes 
considerably smaller than the length-scale for confinement. This division, however, seems 
arbitrary. Indeed, it is not at all clear how and at what density the "transition" should 
be made from a hadronic- to a quark-matter description. Moreover, there might be certain 
regions where none of the pictures could be physically realized. For example, one could 
imagine the existence of an intermediate-density phase characterized by an increase in the 
length-scale for confinement ("swollen" baryons) but also in which clustering correlations 
remain important. Perhaps the most serious difficulty encountered in studying the density 
dependence of strange matter is how to model a system that has quarks confined inside 
color-neutral hadrons at low density but free quarks at high density. The divorce of the 
two pictures is further reinforced by the difficulty of treating quark confinement: how can 
quarks be confined inside hadrons, yet the hadrons separate without generating long-range 
van der Waals forces? We offer no new insight into this difficult problem; rather, we argue 
here that it is useful to consider an effective model which interpolates between a hadron- 
and a quark-based description. 

To do so we rely on the string-flip model. One strength of this model is that the evolution 
of the system from a low-density hadronic phase to a high-density quark phase is dynamical 
rather than through the introduction of artificial parameters. Yet the string-flip model 
employed here is simple and its predictions should be taken with caution. Still, we trust that 
some of the qualitative features obtained from our analysis will have a direct correspondence 
to those made in more sophisticated models. The string-flip model used here is an offspring 
of the original one- dimensional model proposed by Lenz and collaborators |7]] and applied 
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for the first time to nuclear matter by Horowitz, Moniz, and Negele ||. Since then the 
model has been generalized to three dimensions and internal degrees of freedom (such as 
color) have been incorporated to study a variety of ground state properties PHT2|j. All these 



models represent simple — yet realistic — many-body generalizations of the non-relativistic 
constituent quark model ||13| . The cornerstone of the string-flip model is the many-body 
nature of the potential which allows quark confinement within hadrons, while enabling the 
hadrons to separate without generating long-range van der Waals forces JH|. 



There are several questions that we will try to answer in this paper. First, how does the 
length-scale for quark confinement evolve as the density of the system increases and how is 
this evolution affected by the presence of strange quarks? Second, how does the many-body 
(confining) potential modifies the transition to strange matter relative to the predictions of 
a Fermi-gas estimate; do clustering correlations remain important at the transition density 
or has the system dissolved into a collection of noninteracting quarks? Finally, what can one 
learn about the dynamic response of the system through a study of two-body correlations 
between quarks. 

Our paper has been organized as follows. In Sec. |I| we introduce the particular version 
of the string-flip model used here to compute a variety of ground-state properties of strange 
matter. Variational Monte-Carlo results are presented in Sec. |TTT] and compared against 
some useful limiting cases. Finally, our conclusions and outlook are offered in Sec. [IV]. 



II. FORMALISM 



We compute the properties of strange matter in a string-flip model that has constituent 
quarks - - interacting through a many-body potential — as the fundamental degrees of 
freedom. The many-body potential used here is explicitly symmetric in all quark coordinates, 
confines quarks within hadrons, and enables the hadrons to separate without generating long- 
range van der Waals forces. In this paper we use the simplest form of this model that can 
describe the basic quark structure of strange matter. In this version quarks are constrained 
to move in one spatial dimension and have only two SU(2)-like intrinsic degrees of freedom: 
color and flavor. The color — a global, rather than a local — quantum number can be 
either red or blue (with red+blue= "white" ) and the flavor either up or strange. Thus, 
"color-neutral hadrons" are represented as two-quark composites coming in three different 
varieties: nucleon (uu), lambda (us), and cascade (ss). 



A. Many-quark Potential 

The Hamiltonian for a system of A quarks moving in one-spatial dimension under the 
influence of the many-body potential V is given by 

N 2 

H = T + V = Y,£r + V{x). (1) 

The first step in constructing the A-body potential V(x) is the pairing of all A-quarks into 
A = A/2 color-neutral clusters, irrespective of quark flavor. This procedure results in a 
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total of A\ possible configurations. For each of these configurations one now introduces a 
flux-tube energy: 




(2) 



where r\ n and b\ are the positions of the red and blue quarks within the iih cluster, and k 
is the strength of the confining potential (which here has been assumed harmonic). Finally, 
the many-body potential is obtained by demanding that the pairing of quarks into color- 
neutral hadrons be optimal. That is, that the potential be the minimum value of V n among 
all the A\ possible pairings, 



Note that in what follows we adopt units in which h, the spring constant k, and the up-quark 
mass are all equal to one while for the constituent strange-quark mass the canonical value 
of m s /m u = 1.6 is used. 

The minimum prescription described above ensures that the potential is, indeed, sym- 
metric under the exchange of identical quarks. Moreover, no long-range van der Waals forces 
are generated in the model as the force saturates within each individual hadron. In this way 
residual interactions between hadrons involve the possibility of quark-exchange (when the 
strings flip) and the Pauli exclusion principle between identical quarks. Finally, as it is 
possible for all pairings to change as one moves a single quark, the potential is truly an N- 
body operator; it cannot be reduced to a conventional sum of two-body terms. The global 
minimization procedure should be taken among all possible A\ configurations, which rapidly 
becomes unmanageable; for A = 40 (the canonical value used in our simulations) the total 
number of possible configurations is close to 10 48 . Fortunately, because we have restricted 
the problem to one spatial dimension, finding the optimal solution scales linearly — not 
exponentially — with the number of quarks. (In reality the solution to the pairing problem 
scales only with the cube of the number of quarks |ll|Jl5fl . To our knowledge a polynomial 
solution to the three-quark assignment problem has yet to be found). 

To generate the desired A configurations one starts by ordering the red and blue quarks, 
independently, in a one-dimensional box of length L (note that the concept of "order" is 
meaningful only in one spatial dimension). The first of these configurations is obtained by 
mapping the first red quark into the first blue quark, the second red quark into the second 
blue quark, and so on until finally the last red quark is paired with the last blue quark. 
Similarly, the second of these configurations is generated by mapping the first red quark 
into the second blue quark and then continuing the assignments sequentially. One continues 
this procedure until in the last configuration the first red quark is mapped into the last blue 
quark. That is, 




(3) 
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rA^b A , 



1 A 

or V 1 {x) = -Y J {r i 




n <-> b 2 



r 2 «-> b 3 



ta <-> b 



1 A 

or V 2 (x) = ~J2(ri 




i=l 
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n <-> b A ; r 2 <-> 61 ; . . . ; r A <-> 6 A -i , or V^(x) 

The potential energy is then selected to be the minimum among the A configurations: 

V(x) = min {^(x), ^(x), . . . , Vk(x)} . (4) 

Despite the relative simplicity of the model a clear picture of the quark dynamics emerges 
in the low- and high-density phases of the system. In the low-density phase the confine- 
ment scale — a typical distance between quarks within a hadron in free space — is much 
smaller than the inter-hadron separation. Since the potential saturates within each in- 
dividual hadron, quark exchange is suppressed and the system resembles a collection of 
noninteracting nucleons. In the high-density phase, however, the confinement scale becomes 
larger than the inter-hadron separation making string rearrangement dominant. By then 
the potential energy is unimportant relative to the kinetic energy and the system resembles 
a free Fermi gas of quarks. 
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B. Fermi Gas of Quarks 

In this section we compute ground-state properties of a simple Fermi-gas model of quarks 
in order to establish a baseline against which the predictions of the string-flip model may 
be compared. The total energy of a one-dimensional free Fermi gas of quarks is simply the 
sum of the kinetic energies of the two flavors 

^(p,a) = (l-a) + ^(l-ar + rn s a + ^- k F = ^p, (5) 

where the strangeness-per-quark parameter a = \S\/N and the density p have been intro- 
duced. However our results will be presented in terms of the more conventional strangeness- 
per-baryon ratio f s = \S\/A and the relative density p/p c where p c is defined below. Note 
that k F is the Fermi momentum of the system in the limit of |<S| = 0. The strangeness-per- 
quark ratio is determined by demanding that the total energy of the system be minimized 
with respect to a at fixed baryon density. That is, 



'^e)=o = -i-S (1 _ ff) . + m . + ^v. (6) 



s 



The above condition represents chemical equilibrium meaning that the chemical potential 
of both species must be equal. The solution to this quadratic equation is simple and from 
it one obtains the strangeness-per-quark ratio as a function of the density of the system 



m s k F - \\m 



2(m s - l) 2 + k 



a{p) = i^w F ■ (7) 

Although not explicitly stated, the above equation is valid only above a critical density p c 
below which the strangeness-per-quark ratio vanishes. This is due to the fact that below p c 
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it is not yet energetically favorable to introduce strange quarks into the system. The critical 
density is given by 

2 



Pc = -J2{m s -l) ^ 0.697. (8) 

7T v 

Finally, the Fermi-gas energy is obtained by substituting Eq. (0) into Eq. @, i.e., 

T FG (p)=T FG Ua(p)) . (9) 



Another observable useful in characterizing the transition from hadronic- to quark-matter 
is the two-body correlation function. This function represents the ground-state probability 
of finding a pair of quarks at a fixed separation. That is, 

P2(r) = (Y,5(r-\r l -r J \)) = Y,e^ £ (<?U, S AA^) , (10) 

\i<j I 1 ki,si;k2,S2 

where C\ s {Ck, s ) represents a quark creation(annihilation) operator with momentum k and 
(SU(2)-like) color quantum number s. The above two-body operator is easily evaluated in 
the Fermi-gas limit to obtain 

»<0« S «rfr) = l-i(-^)\ (11) 

Note that we have introduced the function gFG{ r ) which is simply the two-body correlation 
function normalized to one at large distances. This function is sensitive only to Pauli 
correlations. The second term in Eq. ([11]) - - with its corresponding 50% reduction at 
the origin — accounts for this fact. 



C. Variational Wavefunction 

In this section we compute ground-state properties of the string-flip Hamiltonian of 
Eq. (|l|) using a variational Monte Carlo approach. For a system of N quarks moving in a 
one-dimensional box of length L the one parameter variational wavefunction has the form 

A = e - xu ^ FG (x) , (12) 

where A is the variational parameter and ^ fg represents a free Fermi-gas wavefunction 
comprised of the product of four Slater determinants; one for each color-flavor combination. 
For example, the one- dimensional Slater determinant for the "red-up" combination is of the 
form IHIBI 



^r UP V) = Ilsm j{r7-rf) . (13) 
i<j lL J 

The Fermi-gas wavefunction is exact for a system of fermions with no correlations other than 
those generated by the Pauli exclusion principle, while the exponential portion of the vari- 
ational wavefunction — symmetric under the exchange of identical quarks — characterizes 
the amount of clustering in the ground state. Here we have introduced the "quasi-potential" 
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U(x) = a uu V uu (x) + a us V us (x) + a ss V ss (x) , (14) 

where V uu , V us , and V ss are the portions of the many-body potential that reside in nucleons, 
lambdas, and cascades respectively; that is: V(x) = V uu (x) + V us (x) + V ss (x). The constants 
Oi U u, ®us, and a ss correct for the sizes of the different hadrons and are thus proportional to 
the ratio of the reduced mass of the corresponding hadron to that of the nucleon: 



a uu = 1 , a us = J ps 1.109 , and a ss = Jm~ s ps 1.265 . (15) 

V 1 + m s 

In the low-density phase — or in a phase with various noninteracting hadrons — the varia- 
tional wavefunction becomes exact in the limit of A — > 1/ \/2 (value of the oscillator param- 
eter for isolated two-quark clusters). Similarly, in the high-density limit where the potential 
is unimportant, the variational wavefunction reproduces exactly that of a free Fermi gas of 
quarks in the limit of A — > 0. 



D. Metropolis Monte Carlo 

An additional advantage of using such a simple variational wavefunction is that the 
expectation value of the kinetic and potential energies are not independent PJIBfl. Indeed, 



the expectation value of the kinetic energy, difficult to simulate because of its derivatives, is 
easily related to the expectation value of the potential energy. That is, 

(* A |T|tt A ) = T FG + 2A 2 (^ A |K|^ A ) , (16) 

where Tfg is the kinetic energy of a one-dimensional free Fermi gas given in Eq. (|5]). The 
minimum value of the kinetic energy is attained in the Fermi-gas limit. The additional term 
in the above expression — 2\ 2 (V)\ — represents the increase in the kinetic energy above 
the Fermi-gas limit due to the presence of clustering correlations. While the kinetic energy 
disfavors such correlations the potential energy favors them, at least at low density. The 
dynamic competition between the kinetic and potential energies will yield an optimal value 
for the variational parameter for any given density and strangeness-to-baryon ratio. The 
optimal value of the variational parameter is thus obtained by minimizing the total energy 
of the system: 

(E) Xpa = T FG (p,a) + (2A 2 + l)(V) Xpa . (17) 

As promised, to compute the ground-state expectation value of the energy we need only 
compute the expectation value of the potential energy, which we accomplish using Monte- 
Carlo methods. Using the algorithm of Metropolis et ai, it can be shown |l7j that the 



expectation value of the potential energy can be obtained by performing a suitable average 
of the form 

i M 

W-^EVM, (18) 

m=l 

where the sequence of points {x m } are distributed according to the square of the variational 
wavefunction. It then becomes a simple — yet computationally-intensive - task to extract 
the optimal variational parameter that minimizes the energy of the system for each value of 
the density and for each number of strange quarks. 
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III. RESULTS 



We start this section by displaying in Fig. [I] the energy-per-quark of the system as a 
function of the variational parameter for a variety of strangeness-per-baryon ratios at a 
fixed density of p = 1.2 (or p/ p c = 1.72; recall that p c = 0.697 represents the critical density 
for the transition to strange matter in a Fermi-gas model). Note that all the calculations 
reported in this work have been effected at the fixed baryon number of A = 40. The plot 
shows the, time-consuming, procedure used to extract the optimal variational parameter, 
the strangeness-per-baryon ratio, and the energy-per-quark at a fixed value of the density. 
Indeed, for every value of f s the optimal variational parameter is extracted as the point at 
which the derivative of the energy with respect to A vanishes. For example, Fig. [I] shows 
that, for this density, the minimum energy-per-quark decreases when the number of strange 
quarks goes from 8 (f s = 0.2) to 24 (f s = 0.6) but increases again if more strange quarks 
(f s = 1.0) are added to the system. The optimal value of A at a density of 1.2 must therefore 
be extracted when the system contains 24 quarks. One then repeats this procedure for all 
densities. 

The optimal variational parameter as a function of density is displayed in Fig. 0. Varia- 
tional results for strange (circles) and non-strange nuclear (squares) matter are shown; the 
solid and dashed lines are smooth interpolations to the data. As expected, the variational 
parameter evolves from the isolated-hadron limit of Ao = 1/ v2 at low density towards the 
Fermi-gas limit of A = at high density. It is important to stress that the decrease in A 
with density — with the corresponding increase in the length-scale for confinement — is a 
prediction of the model. At the transition density to strange matter, p = 0.9 or p= 1.3p c , 
both sets of results start to differ, indeed strange matter favors stronger clustering corre- 
lations because of the larger mass of the strange quark. In this way we obtain one of the 
most important results of our simulations: at the transition density the length-scale for con- 
finement (proportional to A -1 / 2 ) has increased by only 10%, relative to its isolated-hadron 
value. At this density — and well into the strange-matter region — clustering correlations 
remain strong and the system is far from becoming a collection of noninteracting, or weakly 
interacting, quarks. 

The equation of state of the system — the energy-per-quark as a function of density — is 
displayed in Fig. ^. The variational results (circles joined by the solid line) evolve from the 
isolated-hadron to the Hartree-Fock limit. The importance of clustering correlations at low 
density is irrefutable. Indeed, a Hartree-Fock estimate of the ground-state energy (triangles 
joined by the dotted line) diverges at low densities as E HF ~ p~ 3 / 2 . Moreover, although clus- 
tering correlations will eventually become unimportant, the energy of the system remains 
about 5% larger than the Fermi-gas result even at the highest shown density. This picture 
also emphasizes that at the transition region, depicted by the point at which the nuclear 
calculation (squares joined by the dashed line) departs from the strange-matter calculation, 
clustering correlations remain strong. Finally, we have included "hypernuclear" results (dia- 
monds joined by dot-dashed lines) obtained by maintaining the variational parameter fixed 
at its free-space value (A = Ao at all densities). In this approximation baryons are treated as 
incompressible quark composites. It can be seen that for the range of densities probed in 
our simulations, baryon "swelling" — the increase in the length scale for confinement with 
density — seems to play an insignificant role, making the hypernuclear limit an excellent 



S 



approximation to the variational results. 

Even more insensitive to baryon swelling is the strangeness-per-baryon ratio displayed 
in Fig. |j. Indeed, aside from a minor delay in the transition to strange matter relative to 
the Fermi-gas estimate, the strangeness-per-baryon ratio seems insensitive to most details 
of the many-body potential. Clustering is the only feature that seems to affect f s . Yet even 
the Hartree-Fock estimate agrees with the variational result for p> 1.8 p c despite the fact 
that at this density the Hartree-Fock result overestimates the variational energy by almost 
20%. 

We continue by presenting results for the two-body correlation function. While providing 
a useful characterization of the transition from hadronic- to quark-matter, this observable 
also appears to be more sensitive to the fine details of the model. In Fig. ^| we display 
the evolution of the "up-up" correlation function with density. The upper-left-hand panel 
shows the system at very low density where clustering correlations remain strong. Thus, the 
two-body correlation function at short-to-intermediate distances is accurately represented 
by the Gaussian behavior of the wavefunction of an isolated nucleon (depicted by the dashed 
line and barely visible in the figure). Also shown is the two-body correlation function for 
a free Fermi-gas of quarks [solid line; see also Eq. ([11])]. Evidently, Pauli correlations play 
an insignificant role in this very dilute system. As the density increases to p — 0.6, f s 
remains at zero but nucleon swelling — the increase in the length-scale for confinement — is 
now discernible as indicated by the deviation from the isolated-nucleon result. Moreover, the 
correlation function now displays the typical long-distance oscillations, albeit much stronger, 
of a Fermi-gas system. In the third panel (p = 1) the system displays features related to 
both — hadronic and Fermi-gas — limits: an enhanced correlation at short distances and 
a rapid "healing" towards the Fermi-gas limit at larger distances. Finally, the last panel 
shows that the evolution to the Fermi-gas limit is now complete, in spite of a value for the 
variational parameter which, at 50% of its isolated-hadron value, is still large. 

In addition to the up-up correlation function one can measure two-body correlations for 
the other two flavor combinations as displayed in Fig. |[ The upper two panels show, as a 
result of the Pauli principle not being operative between quarks of different flavors, a feature- 
less "up-strange" correlation function. In contrast, the strange-strange correlation function 
at intermediate density (lower-left panel) shows characteristics that cannot be associated 
exclusively with either those of isolated clusters or a Fermi gas but contains some features 
of both. At the largest density shown in the figure (p = 2) the transition to the Fermi-gas 
limit seems now complete. 

So why the emphasis on the two-body correlation function? As stated earlier, the 
two-body correlation function is a useful observable for characterizing the transition from 
hadronic to quark matter. Moreover, it displays clearly the dynamic interplay between clus- 
tering and Pauli correlations. Yet its real appeal stems from its relation to the dynamic 
response of the system. Indeed, the two-body correlation function is related to the Fourier 
transform of the static structure factor: 

—S(q) = 1 +p [°° e~ iqr (g{r) - l)dr = 1 + p f°° 2 cos(gr) (g(r) - l)dr , (19) 

N J-OD v ' Jo v ' 

with the static structure factor, itself, being defined as the integrated sum of the dynamic 
response function, 



9 



S(q) = / S(q,uj)duj. (20) 
Jo 

As a preview of future work, devoted to the study of the exact Euclidean response of the 
system and its analytic continuation into real time, we focus on the static structure factor 
known as the Coulomb sum. In Fig. [7| we display the Coulomb sum as a function of the 
momentum transfer at the low density of p = 0.2. The dashed line represents the response of 
free Fermi gas of quarks; evidently well below the variational result. This behavior suggests 
that at these momentum transfers both quarks in a hadron respond coherently. Yet the 
coherence is incomplete, as the response of the system at low q is below that of a free Fermi 
gas of nucleons (solid line). Some of the coherence is lost because of the intrinsic quark 
structure of the nucleon: the nucleon "survival" probability is proportional to its elastic form 
factor. Note that the very large "spike" at exactly q = 2kp is an artifact of the "rigidity" 
imposed by the Pauli exclusion principle in one spatial dimension. We are currently in the 
process of applying familiar nuclear techniques, such as the impulse approximation, to gain 
insight into the physics of this fundamental problem. Preliminary results for a simplified 
version of the model have been published in Ref. [[18] . 



IV. CONCLUSIONS AND OUTLOOK 

We have employed variational Monte Carlo techniques to compute ground-state prop- 
erties of strange matter in a one-dimensional string-flip model with two flavors and two 
colors of quarks. In this simplified model baryons are color-neutral, two-quark composites. 
We have used a many-body potential that confines quarks within baryons yet enables the 
hadrons to separate without generating unphysical long-range forces as the confining force 
saturates within each individual hadron. The many-body nature of the potential entails a 
difficult optimization problem: how to pair N quarks into A — N/2 color-neutral clusters so 
that the overall flux-tube energy is minimized. For the general three-dimensional case this 
problem involves searching among A\ distinct configurations, a task that becomes prohib- 
itely costly even for small systems. The one- dimensional case studied here is special as one 
can restrict the search to only A configurations. In this manner we were able to simulate 
systems with up to N = 80 quarks (note that in an exhaustive routine one would have to 
search among close to 10 48 configurations!). In this paper we have used, for the sake of 
computational expediency, the simplest string-flip model that can describe the basic quark 
structure of strange matter. 

The crucial feature of all string-flip models is the need to determine an optimal grouping 
of quarks into hadrons. This need alone yields two natural limits for the model: a hadronic 
limit at low density and a Fermi-gas limit at high density. The system thus exhibits a 
transition from hadronic matter, where clustering correlations dominate, to quark matter, 
where only Pauli correlations remain. This is one of the greatest advantage of the model: 
the transition from hadronic- to quark-matter is dynamic, as there is no need to resort to 
ad-hoc prescriptions. The main goal of this paper was the characterization of the transition 
to strange matter; does the transition happens at a density at which the system has already 
dissolved into a collection of "free" quarks or do clustering correlations remain strong? In 
our model the transition to strange matter occurred at a density of approximately p=1.3p c 
(p c was defined as the transition density predicted by the free Fermi-gas estimate). Thus, 
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the many-body potential delays the transition relative to that of a free Fermi gas. More 
interestingly, the length-scale for confinement increased at the transition by merely 10% of 
its isolated-hadron value and in fact clustering correlations remained strong well into the 
transition region. Indeed, we observed that ground-state properties computed in the "hy- 
pernuclear" limit were almost identical to those predicted by the variational approach. In 
sharp contrast, a Hartree-Fock estimate — which includes Pauli but not clustering correla- 
tions — predicted a much larger value for the transition density and an energy-per-quark 
that overestimates (at the transition density) the variational estimate by about 20%. Thus 

- on the basis of our findings — we conclude that the behavior of the system at the tran- 
sition density is well reproduced by a hadronic model with no baryon swelling, but not by 
one that ignores clustering correlations. 

In the future we plan to extend this calculation in several directions. First, we need 
to refine the many-body potential in order to study the possibility of stable, or absolutely 
stable, strange matter. The string-flip potential used here seems to account for the short- and 
intermediate-range structure of the nucleon-nucleon potential. Indeed, the Pauli principle 
between quarks generates a short-range repulsion between clusters while quark-exchange 
generates, at least part of, the intermediate-range attraction. What is missing from the 
model is the long-range (pion) tail. One could introduce a phenomenological long-range part 
to the potential constrained to reproduce the equation of state of nuclear matter. With this 
model for the potential in hand one could then predict the stability, or lack-thereof, of strange 
matter. Second, while undoubtedly challenging, simulating three-dimensional strange matter 
with three flavors and three color of quarks is now within computational reach. Finally, one 
should attempt to calculate the dynamic response of the system. We wish to understand 

- as a function of density and of the momentum transfer — when do leptons scatter from 
hadrons and when do they scatter from the individual quarks. We have offered a preliminary 
answer to this question. By computing the two-body correlation function, and from it the 
Coulomb sum, we observed that the response of the system at small momentum transfers was 
well above the response of a free Fermi gas of quarks. This suggested that all the quarks in the 
hadron respond coherently. Yet we saw that the coherence was incomplete, as the response 
was below that of a Fermi gas of nucleons. We attributed this loss of coherence to the internal 
quark structure of the hadron. How this picture can be reconciled with more sophisticated 
models of the reaction, such as in the impulse approximation, remains an interesting open 
question. In summary, we have shown that in spite of their deceiving simplicity, string-flip 
models of hadronic matter display rich behavior that could yield valuable insights into the 
physics of strange matter. 
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FIGURES 




FIG. 1. Energy-per-quark as a function of the variational parameter for a variety of 
strangenes-per-baryon ratios. The number of baryons has been fixed at A = 40 and the quark 
density at p = 1.2 (p/p c = 1.72). 
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FIG. 3. Energy-per-quark as a function of density. See the text for details. 
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FIG. 5. Two-body up-up correlation function for a variety of densities. 
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FIG. 6. Two-body up-strange and strange-strange correlation functions for a variety of densities. 
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FIG. 7. Coulomb sum as a function of density in units of the Fermi momentum. 
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